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Abstract 

Exclusive semileptonic B decays to the lightest excited charmed mesons are in- 
vestigated at order Aqcjj/mQ in the heavy quark effective theory. At zero recoil, 
Aqcd/™q corrections to the matrix elements of the weak currents can be writ- 
ten in terms of the leading Isgur-Wise functions for the corresponding transition 
and meson mass splittings. The differential decay rates are predicted, including 
Aqcd/™Q corrections with some model dependence away from zero recoil. Appli- 
cations to B decay sum rules and factorization are presented. 
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1 Introduction 

Heavy quark symmetry S implies that in the hmit where the heavy quark mass 
goes to infinity, matrix elements of the weak currents between a B meson and 
an excited charmed meson vanish at zero recoil (where in the rest frame of the 
B the final state charmed meson is also at rest). However, in some cases at 
order Aqcb/ttt-q these matrix elements are not zero.B Since most of the phase 
space for semileptonic B decay to excited charmed mesons is near zero recoil, 
^QCD /"niQ corrections can be very important. 

In the mg — > oo limit, hadrons containing a single heavy quark, Q, are 
classified not only by their total spin J, but also by Ihc spin of their light 
degrees of freedom (i.e., light quarks and gluons), si.a Such hadrons are in 
degenerate doublets with total spin J± = s; ± i, coming from combining the 
spin of the light degrees of freedom with the spin of the heavy quark, sq — ^. 

The ground state mesons contain light degrees of freedom with spin-parity 
s^' = i , giving a doublet containing a spin zero and spin one meson. For 
Q = c (b) these are the D and D* [B and B*) mesons. Excited charmed 
mesons with ' = | have been observed. These are the Di and D2 with 
spin one and two, respectively. In the nonrelativistic constituent quark model 
these states correspond to L = 1 orbital excitations. Combining the unit of 
orbital angular momentum with the spin of the light antiquark leads to states 
with s^' = i"*" and l"*^. The doublet, {D^,D\), has not been observed, 
presumably because these states are much broader than those with s^' = 
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The hadron masses give important information on some HQET matrix 
elements. The mass formula for a spin symmetry doublet of hadrons H± with 
total spin J± = s; ± i is 

where n± = 2J± + 1 is the number of spin states in the hadron H±. A is the 
energy of the light degrees of freedom in the mg oo limit, Ai determines the 
heavy quark kinetic energy, and A2 determines the chromomagnetic energy. 
These parameters depend on the particular spin symmetry doublet to which 
H± belong. We reserve the notation A, Ai, A2 for the ground state multiplet 
and use A', A'^, X'2 and A*, A*, for the excited sj' = and doublets, 
respectively. The excitation energy A' — A plays a very important role, 

-^,_-^^ m, {m'^ ~rnn)-m. {rn'^ " + f ^\ ^ 0.39 GeV , (2) 
mb- rric \ J 

where Ttih — {n^niH^ + n^mH_)/{n^ + n_) are the spin average masses. 
This value of A' — A has only a small sensitivity to mt and nic (we used 
nib = 4.8 GeV and rric = 1.4 GeV), but it has considerable uncertainty because 
the experimental error on rn'^ is large. (We identified the B^*^7r resonances 
observed at LEP with the bottom s^' — |^ meson doublet, to^ = 5.73 GeV.@) 

2 The B DieVf, and B D2eV(, decays 

The matrix elements of the vector and axial- vector currents {V^ = 07^ 6 and 
= cj^j5 b) between B mesons and Di or D2 mesons are parameterized as 



{D,iv',e)\V'^\Biv)) 



^niDi TUB 
(Z.,(.-,e)|AMi.(.))^^^^^^^ 

^TODi riiB 
{D*2W,e)\A'^\B{v)) 
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{Dl{v\e)\V^^\B{v)) 



^ikve^^''P'<el,v-vpv' (3) 



where the form factors fi and ki are dimensionless functions oi w ~ v ■ v' . At 
zero recoil {v — v') only the fvi form factor can contribute, since v' dotted 
into the polarization (e*'' or e*'^") vanishes. 



To write the differential decay rates in terms of the form factors in Eq. ^ 
we define 9 as the angle between the charged lepton and the charmed meson in 
the rest frame of the virtual W boson, i.e., in the center of momentum frame 
of the lepton pair. The different helicities of the Di or D2 yield distinct dis- 
tributions in 9, which will probably be measurable. These helicity amplitudes 
are affected by the 1/toq corrections in different ways, so measuring them will 
be interesting. In terms oi w — v ■ v' and cos 9, the differential decay rates are 

dTn 



3rl 1 sin2 9 [{w - n)fv, + {w^ - l)ifvs + nfv^) 

+ (1 - 2nw + rl) \{1 + cos^ 9) + {w^ - 

dw dcos 9 



dr 



-Do* 



To 



2rl (u)2 



1 



,3/2 



4:Cos9\/w^ - Ja 



(4) 



- (1 - 2r2W + rl) yil + cos^ 9) [k\^ + {w^ - l)fc^] 



4 cos 9 \/u)2 — 1 kA-^ ky 



dw dcos 9 , 



(5) 



where Tq — Gp\Vcb\'^m^/{192TT^), ri — mD^/niB and r2 = mD^/rriB- The 
semileptonic B decay rate into any J 7^ 1 state involves an extra factor of w^ — 1. 
The sin^ 9 term is the helicity zero rate, while the 1 + cos^ 9 and cos 6* terms 
determine the helicity A = ±1 rates. Since the weak current is ^ — A in the 
standard model, B mesons can only decay to the helicity |A| = 0, 1 components 
of any excited charmed mesons. The decay rate for |A| = 1 vanishes at maximal 
recoil, Wmax = {l+r^)/{2r), as implied by the l — 2rw+r'^ factors above (r = ri 
or r2). The differential decay rates in terms of the electron energy follow from 
Eqs. ^ and ^ using y = 2Ee/mB = 1 — rw — rV — 1 cos 9. 

The form factors fi and ki can be parameterized by a set of Isgur-Wise 
functions at each order in AqcD/^^Q.—In the toq ~* 00 limit, fi and ki are 
given in terms of a single function r(K;).EI Only fvi can contribute to the matrix 
elements in Eq. |^ at zero recoil. Heavy quark symmetry does not fix t(1), since 
/vi (1) = in the infinite mass limit independent of the value of r(l). 

At order 1/toq several new Isgur-Wise functions occur. Corrections from 
matching QCD currents onto HQET currents introduce dependence on A', 
A, and two new functions, ti.2(w). The order 1/toc chromomagnetic correc- 
tion to the Lagrangian introduces three new functions, 771,2,3 (iw), while the 
charm quark kinetic energy term introduces dependence of the form factors 
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on one more function, 77ko(w). The functions parameterizing 1/mt, correc- 
tions to the Lagrangian occur in a single hnear combination, rjb (w) . Denoting 
Ec.fc — ^/^uicb, the B — > Dieve form factors including order l/rric^b terms area 

%/6/a = -{w + 1)t - Ebliw - 1)[(A' + A)t - {2w + l)ri - T2] + {w + l)7jb} 
-£c[4(wA' - A)t - 3{w - 1)(ti - ra) + {w + l){Tjkc - 2m - 3m)] , 
V6fv^ = (1 - w^)t - Sbiw^ - 1)[(A' + A)r - (2w + l)ri - rs + 77^] 

-£,[4(u; + l)(w;A' - A)t - (u;^ - 1)(3ti - 3r2 - 7?ke + 2r/i + 3%)] , 
Vefv, - -3t - 3e6[(A' + A)t - (2i(; + 1)ti - T2 + j/t] 

-£c[(4w - l)ri + 5t2 + 377kc + Wvi + 4(w - 1)772 - 5%] , 
V^/V3 = (w' - 2)t + £fc{(2 + w)[{A' + A)t - (210 + 1)ti - r2] - (2 - w)7^b} 
+ec[4(wA' - A)r + (2 + w)ti + (2 + 3u))t2 

+ (w - 2)77ko - 2(6 + u>)?7i - 4(w - 1)772 - (3w - 2)7/3] • (6) 

The analogous formulae for B Z^j^^^e are 

kv = -T - e6[(A' + A)t - (2u; + 1)ti - T2 + 775] 

-ec(Ti - r2 + 77ko - 27/1 + 7/3) , 
fcAi = -(1 + w)r - Ebiiw - 1)[(A' + A)t - (2u; + 1)ti - T2] + (1 + u;)7yb} 
-ec[(7« - 1)(ti ~ T2) + (w + l)(77kc - 2771 + 773)] , 
= -2ec(Ti + 772) , 
kAs = T + £b[(A' + A)t - (2u; + 1)ti - T2 + r]b] 

-Scin + T2 - 77kc + 2771 - 2772 - 773) . (7) 

At order l/mq the form factor fvii^) is no longer zero. At this order it 
can be written in terms o£ A' — A and the Isgur-Wise function T(7i;) evaluated 
at zero recoil. Explicitly, □ 

V6/y,(l)--— (A'-A)t(I). (8) 

The factor of four in the numerator of Eq. |^ makes this quite a large correc- 
tion. Furthermore, its importance is enhanced over other Aqcd/^^iq correc- 
tions since most of the phase space is near zero recoil. (For a flavor diagonal 
current a relation similar to Eq. ^ was previously obtained by Voloshin. a) 

The allowed kinematic range for B DieD^ [B DjcS'e) decay is 1 < 
w < 1.32 {1 < w < 1.31). Since these ranges are fairly small, and there are 
some constraints on the l/mg corrections at zero recoil, it is useful to consider 
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the decay rates given in Eqs. ^ and |^ expanded in powers of w — 1. The general 
structure of the expansion of dT/dw is elucidated schematically below, 

.(A=0) 



dr 



dw 



+ (u>-l)2(l + e + . 



dr 



dr 



dw 



(|A|=0,1) 
D' 



(0 + 
+ {w 



Oe 



1)2(1 



...) + (u'-l)(0 + £ + £^ + ...) 
,.) + .. 

...) + (u;-l)(l + £ + ...) 
..) + .. 



dw 



[W^ ~ 1)3/2 \^l + s+...) + {w- 1)(1 +£ + ...) + 



(9) 



Here e" denotes a term of order (AqcD/Tig)". The zeros in Eq. ^ are conse- 
quences of heavy quark symmetry, as the leading contribution to the matrix 
elements of the weak currents at zero recoil is of order Aqcd/™q- Thus, the 
Di decay rate at w = 1 starts out at order Aq(-.j-,/toq. Similarly, the vanish- 
ing of fvii^) in the tuq — > oo limit implies that at order w — 1 the Z?J^~°^ 
rate starts out at order Aqcd/™q- The I?2 decay rate is suppressed by an 
additional power of w2 — 1, so there is no further restriction on its structure. 



3 Applications 

3.1 Predictions for B — > {Di, D2)ev^ decays 

Predictions for various quantities can be made assuming a linear dependence 
of the Isgur-Wise function, t{w) = t(1)[1 + t' {w— 1)]. Lagrangian corrections 
from the chromomagnetic term are expected to be small, so 771,2,3 are neglected 
in Eqs. ^ and 0. In these equations 77^0 and rji, are also set to zero, since the 
contributions from the kinetic energy operator can be absorbed into t. Then 
the unknown parameters entering our predictions are f', ti and T2. Approxi- 
mations Bi and B2 correspond to two choices of ti and T2 (ri = T2 = for Bi; 
Ti = At and T2 — —A' r for B2). The difference between these approximations 
gives a rough estimate of the uncertainty due to unknown 1/toq corrections. 

Recently the ALEPHQ and CLEOB Collaborations measured, with some 
assumptions, the B — > DieDe branching ratio. The average of their results is 

^(B^DiePe) = (6.0 ±1.1) X 10-3. (10) 

The B D2 e De branching ratio has not yet been measured; CLEO and 
ALEPH set upper bounds on B{B D2ei>e) at the 1% level. 
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Table 1: Predictions for the ratio of B — > D\ev^ and B DJefe decay rates, and the 
extracted value of t(1). These results correspond to f' = T'(l)/r(l) = —1.5. 



Approximation 


R = Td'JTd, r(l) 


6.0 X 10-3 


1/2 


B{B -> Dieve) 


mq — !■ oo 
Finite niQ | 


1.65 1.24 
0.52 0.71 
0.67 0.75 



Tablcj^ presents R — T d*/T Bi and the value of t(1) extracted from Eq. |l0|, 
both in the infinite mass hmit where Sc = £b = ^ and for the physical values 
of Ec = 1/(2.8 GeV) and = l/(a^pi|lc)n. In this table f'(l) = -1.5 is used, 
motivated by model calculations. ltE^IOO The 1/toq corrections enhance the 
Di rate significantly, while they affect the D2 rate to a much smaller extent. 
Thus, the values of r(l) and R are substantially reduced from their values in 
the infinite mass limit. R is fairly insensitive to f' and T2, but depends strongly 
on Ti. r(l) is insensitive to T2, but f' and ri affect it at the 20% level. For a 
more detailed discussion and other predictions, see Ref. [2]. 

With more experimental information on the differential decay rates for 
B — > DieDe and B D2ei>e it should be possible to determine from experi- 
ment the values of f ' and ri . This would remove much of the uncertainty in 
the predictions presented in Table |l|. 

3.2 Factorization 

Factorization should be a good approximation for B decay to charmed mesons 
and a charged pion. Neglecting the pion mass, the two-body decay rate, F^, is 
related to the differential decay rate dFsi/dw at maximal recoil for the semilep- 
tonic decay with the tt replaced by the ePg pair. This relation is independent 
of the identity of the charmed meson in the final state, 

1 TT — 2 \ ~] 

r \ aw 

Here r is the mass of the charmed meson divided by uib, Wmax = (1 + ''^)/(2r), 
and /tt ~ 132 MeV is the pion decay constant. C is a combination of Wilson 
coefficients of four-quark operators, and numerically ClKidl — 1- These non- 
leptonic decay rates can therefore be predicted from dFsi/diy at maximal recoil. 
In the absence of a measurement of the differential decay rates, we can use our 
results for the shape of dFgi/dw to predict the B Ditt and B — > rates. 
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At present there are only crude measurements of the B{B Ditt) and 
B{B D^tt) branching ratios. Assuming B{Di{2420f^ D*+tt-) = 2/3 and 
B{D*{2Amf -> D*+TT-) = 0.2, the measured rates areB 

^ L»i(2420)°7r-) = (1.17 ± 0.29) x 10-^ 
B{B- i:>*(2460)°7r-) = (2.1 ± 0.9) x lO'^ . (12) 

Our prediction for B{B — > Ditt)/B{B DieVf.) varies between 0.13 and 
0.31 depending on f', fairly independent of ri^2- Assuming that factoriza- 
tion works at the 10% level, a precise measurement of the B{B -Divr) 
rate may provide a determination of f'. The present experimental data, 
B{B Ditt)/B{B — > Dieve) — 0.2, does in fact support f' '--^ — 1.5, which we 
took as the "central value", motivated by model calculations. □ 

The prediction for B{B D2Tt)/B{B Din), on the other hand, only 
weakly depends on f', but it is more sensitive to T2 and especially ti. Varying 
^1,2, we can accommodate almost any value of B{B D2Tt)/B{B Din) 
between and 1.5. In the future, experimental data on this ratio and R may 
lead to a determination of fi. 



3.3 B {D^,Dl)ePe decays 

Similar results to those presented in Sec. 2 also hold for decays to the sj' = i"*" 
doublet, B (Dq, D^)ez/e.H The zero recoil matrix elements are determined 
by the excitation energy. A* — A ~ 0.35 GeV, and the toq oo Isgur-Wise 
function for these transitions, ({w), at zero recoil. Numerically, the D^ rate 
is enhanced similar to the Di rate, whereas l/mg corrections to the Z?* rate 
enter near w = 1 proportional to the anomalously small combination Sc — 3£6. 

To obtain even a crude absolute prediction for the B D\,Dq rates, a 
relation between the s^' = ^ and | Isgur-Wise functions is needed. In any 
nonrelativistic constituent quark model with spin-orbit independent potential, 
C and T are related byOcI 

C{w) = ^T{w), (13) 

since both the ^ and | doublets correspond to L = 1 orbital excitations. 

This relation implies that the six lightest charmed mesons contribute ahmit 
8.2% of the B decay rate [D and D* is about 6.6% of the total B decay ratellj). 
Therefore, semileptonic decays into higher excited states and non-resonant 
multi-body channels should be at least 2% of the B decay rate, and possibly 
around 3% if the semileptonic B branching ratio is closer to the LEP result of 
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about 11.5%. Such a sizable contribution from higher mass charmed mesons 
and non-resonant modes would soften the lepton spectrum, and may make 
the agreement with data on the inclusive lepton spectrum worse. If ^ were 
enhanced by a factor of two compared to Eq. |l3|, then semileptonic B decay 
rate to the six lightest charmed mesons could add up to close to 10%. However, 
model calculationsc£l seem to obtain a suppression rather than an enhancement 
of C compared to Eq. |l^. Thus, taking the measurements for the B ^ D, D*, 
and Di semileptonic branching ratios on face value, a decomposition of the 
semileptonic rate as a sum of exclusive channels seems problematic. 



3.4 Sum rules 

Our results are important for sum rules that relate inclusive B -^<t^e decay 
to the sum of exclusive channels. For example, the Bjorken sum rulell^ bounds 
the slope of the B D^*^ei'e Isgur-Wise function, defined by the expansion 
^(w) — 1 — {w ~ 1) + . . . . Knowing would reduce the uncertainty in the 
determination of \Vcb\ from the extcapolation of the B — > D^*^eDe decay rate 
to zero recoil. This sum rule reads E3i3 

m p 

The ellipses denote contributions from non-resonant channels. C^™^ and r'^' 
are the Isgur-Wise functions for the exited sf' = and I"*" states, respec- 
tively. (For m = p = these are the lightest orbitally excited states, and 
m,p > I are radial excitations of these.) Since all contributions to the right- 
hand-side of Eq. |l^ are non-negative, a lower bound on can be obtained by 

keeping only the first few terms. Using t(1) ~ 0.71 and Eq. |l3|, we find a large 

1 + T + 

contribution to the Bjorken sum rule from the lightest ^ and ^ doublets, 

,^>i + fil)!! + 2Ml)!!=,„.„. (15) 

4 4 3 

The contribution of the states to this bound, which relies on the quark 
model result in Eq. |l^, is only 0.17. 

The contribution of excited states to other sum rules is discussed in Ref. [2] . 



4 Conclusions 



• At zero recoil, order 1/toq contributions to semileptonic B Di,D2 
decays (any excited charmed meson with + parity) are determined by 
the mg 00 Isgur-Wise function and known hadron mass splittings. 
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• Shape of decay spectra predicted near zero recoil, including l/mg cor- 
rections with reasonable assumptions. 

• Large 1 /mg corrections to some predictions can be tested against exper- 
imental data in the future. 

• Better understanding of inclusive = ^ exclusive in semileptonic B decay. 

• Test applicability of heavy quark symmetry for B decays to excited 
charmed mesons. 
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